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Neutral beam model for the anomalous 7-ray emission 
component in GRB 941017 
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Abstract. Gonzalez et al. (2003) have reported the discovery of an anomalous radiation component from ~ 1 - 
200 MeV in GRB 941017. This component varies independently of and contains xL 3x the energy found in the 
prompt ~ 50 keV - 1 MeV radiation component that is well described by the relativistic synchrotron-shock model. 
Acceleration of hadrons to very high energies can give rise to two additional emission components, one produced 
inside the GRB blast wave and one associated with an escaping beam of ultra-high energy (UHE; Si 10 14 eV) 
neutrons, 7 rays, and neutrinos. The first component extending to ~ 100 MeV is from a pair-photon cascade 
induced by photomeson processes with the internal synchrotron photons coincident with the prompt radiation. 
The outflowing UHE neutral beam can undergo further interactions with external photons from the backscattered 
photon field to produce a beam of hyper-relativistic electrons that lose most of their energy during a fraction of a 
gyroperiod in the assumed Gauss-strength magnetic fields of the circumburst medium. The synchrotron radiation 
of these electrons has a spectrum with vF v index equal to +1 that can explain the anomalous component in 
GRB 941017. This interpretation of the spectrum of GRB 941017 requires a high baryon load of the accelerated 
particles in GRB blast waves. It implies that most of the radiation associated with the anomalous component is 
released at ^ 500 MeV, suitable for observations with GLAST, and with a comparable energy fluence in ~ 100 
TeV neutrinos that could be detected with a km-scale neutrino telescope like IceCube. 
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1. Introduction 

Based on joint analysis of BATSE LAD (Large Area 
Detector) and the EGRET TASC (Total Absorption 
Shower Counter) data, Gonzalez et al. (2003) recently re- 
ported the detection of an anomalous MeV emission com- 
ponent in the spectrum of GRB 941017 that decays more 
slowly than the prompt emission detected with the LAD in 
the « 50 keV - 1 MeV range. The multi-MeV component 
lasts for 200 seconds (the tgo duration of the lower- 
energy prompt component is 77 sec) , and is detected with 
the BATSE LAD near 1 MeV and with the EGRET TASC 
between « 1 and 200 MeV. The spectrum is very hard, 
with a photon number flux <p(e) oc e , where e = hv /m e c 2 
is the observed dimensionless photon energy. 

This component is not predicted or easily ex- 
plained within the standard leptonic model for GRB 
blast waves, though it possibly could be related to 
Comptonization of reverse-shock emission by the for- 



ing self-absorbed reverse-shock optical synchrotron radi- 
ation (Pe'er and Waxman (2004)|). Another possibility is 



ward shock electrons < |Granot and Guetta (2003) I, includ- 
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that hadronic acceleration in GRB blast waves could be 
responsible for this component. 

We propose a model involving acceleration of hadrons 
at the relativistic shocks of GRBs. A pair-photon cas- 
cade initiated by photohadronic processes between high- 
energy hadrons accelerated in the GRB blast wave and 
the internal synchrotron radiation held produces an emis- 
sion component that appears during the prompt phase. 
Photomeson interactions in the relativistic blast wave also 
produce a beam of UHE neutrons, as proposed for blazar 



jets (Atoyan and Dermer (2003) I. Subsequent photopion 
production of these neutrons with photons outside the 
blast wave produces a directed hyper-relativistic electron- 
positron beam in the process of charged pion decay and 
the conversion of high-energy photons from 7r° decay. 
These energetic leptons produce a synchrotron spectrum 
in the radiation reaction-limited regime extending to ^ 
GeV energies, with properties in the 1 - 200 MeV range 
similar to that measured from GRB 941017. If our model 
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Fig. 1. Photon energy fluence from an electromagnetic 
cascade initiated by photopion secondaries in a model 
GRB, with parameters given in the text. Five generations 
of Compton (heavy curves) and synchrotron (light curves) 
are shown. The first through fifth generations are given by 
solid, dashed, dot-dashed, dot-triple-dashed, and dotted 
curves, respectively. The total cascade radiation spectrum 
is given by the upper bold dotted curve. 

is correct, detection of this component therefore gives im- 
portant evidence for the acceleration of UHE cosmic rays 
in GRB blast waves, and would be confirmed by the de- 
tection of UHE neutrinos from GRBs. 

2. Hadronic Model for GRB 941017 

We assume that efficient proton acceleration to ultra- 
relativistic energies takes place in GRB blast waves, as re- 
quired in models where GRBs accelerate high-energy cos - 
mic rays ( |Vietri (1995)1 |Waxman (1995)| permer (2002) I . 
Photopion interactions of these protons with internal syn- 
chrotron photons and with photons from an external radi- 
ation field create neutral particles and charged pions that 
decay and initiate an electromagnetic cascade within the 
GRB blast wave. 

During the prompt emission phase, an additional com- 
ponent of radiation is produced by the cascade occur- 
ring within the blast wave. Fig. 1 shows the hadron- 
initiated cascade radiation for a model GRB at red- 
shift z = 1, with hard X-ray fluence $tot — 3 x 
10 _4 erg cm -2 , a light curve of 100 second dura- 
tion divided into 50 pulses of 1 second each, and 
Doppler factor 5 = 100 (see Atoyan and Dermer (2003) 



Dermer and Atoyan (2003) for more details about the 
model). The total amount of accelerated proton energy 
E' = 47rd^$ tot (5~ 3 (l + z) _1 is injected into the comoving 
frame of the GRB blast wave. 

The various generations of synchrotron and Compton 
radiation initiated by the cascade are shown in Fig. 1, 
along with the total radiation spectrum. As can be seen, 
the isotropic cascade radiation approaches the spectrum 



of an electron distribution cooling by synchrotron losses, 
that is, a photon number index between —1.5 and —2. 
During intervals (a) and (b) of GRB 941017 when the 
prompt hard X-ray and soft 7-ray flux is brightest and 
the high-energy spectral index is hard though not well- 
measured (Gonzalez et al. (2003) I, this component could 
make the dominant contribution to £ 1 MeV radiation 
fluence. The isotropic cascade radiation will decay at the 
same rate as the synchrotron radiation, after which a sec- 
ond component from the ouflowing neutral beam begins 
to dominate. 

Ultra-rclativistic neutrons formed in the reaction p + 
7 — > n + 7T° arc not confined by the magnetic field in 
the GRB blastwave shell and flow out to create an ener- 
getic neutron beam. These neutrons are subject to fur- 
ther photopion processes with photons in the surrounding 
medium to form charged and neutral pions. In the Gauss- 
type magnetic fields surrounding GRB sources that we 
assume here, charged 7r and fx at energies ^ 10 18 eV de- 
cay rather than lose energy through synchrotron emission 
(Rachen and Meszaros (1998) I. The charged pions decay 
into ultrarelativistic electrons and neutrinos, whereas the 
decay of ir° produces two 7 rays that are promptly con- 
verted into electron-positron pairs on this same exter- 
nal radiation field. These energetic electrons (including 
positrons) are initially produced in the direction of the 
GRB jet. 

The spectra of secondary electrons created by the neu- 
tron beam displays a sharp cutoff at energies £ 10 14 eV as 
a consequence of the high threshold for photomeson inter- 
actions (see calculations in Dermer and Atoyan (2003) I. 
Electrons with Lorentz factor 7 lose energy through syn- 
chrotron radiation in an ordered magnetic field with 
strength B at the rate —dj/dt = ot-B 2 7 2 sin 2 ip / (47rm e c) , 
where ip is the electron pitch angle. The corresponding 
synchrotron energy- loss time scale t syn — j/\dj/dt\. The 
gyration frequency lob = eB/jm e c, and is independent of 
pitch angle. When uj B t S y n <c 1, the electron loses almost 
all of its energy into synchrotron radiation in a time less 
than the gyroperiod. We use the term "hyper-relativistic" 
to refer to electrons in this radiation-reaction regime 
of synchrotron emission ( Nelson and Wasserman (1991)[ 
Aloisio and Blasi (2002)1. 



Electrons which cool before being deflected by an an- 
gle 9 equal to the jet opening angle 6j will emit most of 
their energy within 9j . The pitch angle if) does not change 
in the process of synchrotron losses when 7 > 1. Then 
cos 9 = cos 2 if) + sin 2 -0 cos (f>, where </> = uist is the ro- 
tation angle. In the limit of small 9 and cf>, 9 = fismip. 
The condition 9 < 9j for times t < t syn then results 
in the condition 7 ^ 7hr(^j) = \J Amej(9j<j^B sin-0) — 
3 x 10 8 /[ V /(6» J - /0.1)S(G)sin^] for the hyper-relativistic 
electrons. Here we have taken a typical jet opening half- 
angle 9j =0.1 because this is the average value implied by 
the analysis leading to the standard energy reservoir re- 
sult of GRBs ( |Frail et al. (2001)1 . Lower-energy electrons 
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with 7 < jhr(Qj) radiate their energy over a much larger 
solid angle. 

The mean energy E 1 = m e c 2 t — Ej of synchrotron 
photons emitted by electrons with 7 ~ Jhr(9j) is indepen- 
dent of tp and B, and is given by 



The synchrotron emission spectrum at energies E~ -C 



E, 



heB simp -tl r {9j) 



500 



,c (1 + z) (^/0.1)[(l + z)/2] 



MeV, (1) 



Hyper- relativistic electrons with 7 > r jhr{9j) rapidly lose 
energy through synchrotron losses and deposit all of their 
energy along the direction of the jet. Electrons at lower 
energies are deflected to angles 9 > 9j, and their emission 
is not seen by an on-axis observer. Hence the distribution 
of electrons along the jet direction always has an effective 
low-energy cutoff at 7^ r (0j). The production spectrum of 
the electrons can have an intrinsic cutoff - lco clue either to 
the low-energy cutoff in the escaping neutron spectrum, or 
to the neutron-induced photomeson secondary spectrum 
outside the GRB blast wave. If p = j C o/jhr(9j) > 1, 
then the observed synchrotron spectrum is a power law 
with —1.5 index for Ej ^ E 1 ^ p 2 Ej, and a photon spec- 
trum with the same spectral index as the accelerated pro- 
tons and escaping neutrons at photon energies E n ^ p 2 Ej 
1 Nelson and Wasserman (1991) I. If p < 1, then the ob- 
served photon spectrum at E 1 Ej has the same spectral 
index as the primary neutrons. 

At photon energies E 1 <C Ej, the observed spectrum 
is produced by the same hyper-relativistic electrons with 
7 ^ Jhr{9j), but at energies e well below the peak energy 
3j 2 cb, where es = B/B cr and B cr = 4.41 x 10 13 G is the 
critical magnetic field. We now derive this spectrum. 

The differential energy radiated per dimensionless en- 
ergy interval de per differential solid angle element dfl in 
the direction 9 with respect to the direction of an electron 
moving with Lorentz factor 7 is given by 



dE 



dedfl 37r 2 A c jcb 



( ) 2 (1 + 7 2 2 ) 2 (A||+A ± ), 



(2) 



where Ac = Ti/m e c — 3.86 x 10 11 cm is the elec- 
tron Compton wavelength, and Ay = ^2/3 (0 anc ^ ^-L = 
(-f9) 2 Kf/ 3 (£ i )/{l + (j9) 2 ] are factors for radiation polar- 
ized parallel and perpendicular to the projection of the 
magnetic field direction on the plane of the sky defined 
by the observer's direction (Jackson (1975) I. The factor 
£ = e/e, where e = 3e B 7 2 /(l + i'ff'Y'*, and K n (x) is a 
modified Bessel function of the second kind, with asymp- 
totes K n (x) — > \T(n)(2/x) n in the limit x <C 1, and 
K n (x) — > y/ir/2x exp(— x) in the limit x 3> 1. The con- 
dition £ <C 1 corresponds to e <C e where K n (£) are in 
their power-law asymptotes, and £ ^ 1 or e £ e is where 
if n (£) are in exponential decline. The characteristic en- 
ergy e approaches 3es7 2 when 7^ < 1, and e declines 
with 9 according to the relation e = SeB'f 2 / {j9) 3 when 
j9 > 1. When e < e, then A|| > Aj_ and dE/dedVt ~ 
3 1 / 3 (1.07e/7r) 2 (7e/e B ) 2 / 3 /A c oc e 2 / 3 . For a fixed value of 
e, this emissivity exponentially cuts off when e ^ 3es/j9 3 , 
or when 9 > 9 max = (3es/7e) 1/3 . 



Ej, integrated over solid angle, is thus 



dE n r&Bhtf" dE 



dc 







if! 



(3) 



This differs from the energy index +1/3 for synchrotron 
radiation radiated by electrons in the classical regime, be- 
cause in this case one integrates over a complete orbit of 
an electron (Rybicki and Lightman (1979) I. In this case, 
the solid angle element d£l — > <i(9sin# — > d9 simp in the 
integration in eq. J3J). 

3. Neutron Beams from GRBs 

In our model, ultrarelativistic protons undergo pho- 
tomeson interactions with internal and external radia- 
tion photons, producing a beam of outflowing neutrons. 
Subsequent interactions of these neutrons with photons 
of the external radiation field generate a beam of hyper- 
relativistic electrons. The £ 1 MeV emission during the 
prompt phase is mainly cascade radiation produced within 
the GRB shell, while the emission induced by the neu- 
tron beam is formed at later times. This interpretation 
requires that the fluence of photohadronic cascade radi- 
ation within the shell exceed by a factor of order unity 
the cascade fluence induced by the neutral beam radia- 
tion, consistent with observations. At later times, as shown 
above, synchrotron radiation from the hyper-relativistic 
electrons has a low-energy cutoff at _E 7 < Ej ~ 0.5 
GeV, with a specific characteristic power-law number 
spectral index equal to —1, as observed for GRB 941017 



( Gonzalez et al. (2003) I. 

The w 200 s decay time of the anomalous emission can 
be explained by the projected thickness ~ R(l — cos9j) 
of the emitting region from the front to the edges of a jet 
blastwave at distances R w 6 x 10 14 (6» J /0.1)" 2 [(l + z)/2]- 1 
cm from the central source, implying significant opacity 
to photomeson processes due to external photons at these 
distances. Note that hyper-relativistic electrons can make 
a substantial contribution to emission at 9 « 9j ^> 1/F, 
where T is the GRB blastwave Lorentz factor. 

The external radiation field could be due 
to plerionic emission in the supranova model 
(Konigl and Granot (2002) 1, or to GRB synchrotron 
photons that are backscattered by stellar wind electrons 
in the collapsar model ( Beloborodov (2002) I . We estimate 



the effective photomeson energy-loss scattering depth 
t„ 7 for neutrons in the collapsar-model case, and show 
that it implies a size scale of the back-scattering material 
in accord with the radius inferred from the duration of 
the anomalous component. Above the photopion thresh- 
old for backscattered radiation, which applies to 
several hundred GeV protons, r„ 7 w K ni a ni <^ 1 (R)TT, 
where the product of the inelasticity K n ^ and cros s 



section er„ 7 is = 70 ^barns (Atoyan and Dermer (2003)1, 
$ 7 (i?) = d 2 L if/[(e)m e c 2 R 2 } is the photon fluence of the 
prompt GRB emission at radius R, and tt — <Jt:Rii w (R) 
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is the Thomson depth for scattering photons with 
observed energies (e)m e c 2 at radius R = 10 14 i?i 4 cm 
in the frame of the stellar wind. The expression for 
photon fiuence is written in terms of the energy fluence 
if — 10~ 4 </?_4 ergs cm -2 from a GRB at luminosity 
distance d,L = 10 28 d28 cm. 

The stellar wind outflow rate M = 10~ 4 M_ 4 Af Q /yr 
and wind speed v = 2000i>2ooo km s _1 are scaled 
to the properties of the winds of Wolf-Rayet stars 



an IceCube-class km-scale telescope from this type of 
event, in agreement with phcnomcnological inferences by 



Alvarez-Muniz, et al. (2003) 



( Lozinskaya (1992) I, which are likely GRB progenitors. 
This implies a Thomson depth tt = 0.01Af_ 4 /i? 14 ii2oooj 
so that r„ 7 » d| g ( y 9_ 4 A/_4/[(e)(l + z)Ri 4 v 2 ooo] at scales 
R. Equating the restriction on radius following from 
the requirement r„ 7 ^ 1 with the relation on R from 
duration for z = 1 and e?28 = 2.2 implies 9j = 
0.14[((e)/0.1)«2ooo/^4M_ 4 ] 1/6 , which is a typical jet 
opening angle. The differing GRB external radiation and 
density environments which determine the intensity of tar- 
get photons could account for the unusual spectrum of 
GRB 941017. 

The fluence in the anomalous MeV component of GRB 
941017 is £ 3x larger than the prompt emission fluence. 
If the anomalous fluence is assumed to originate from pho- 
topion processes from nonthermal protons, then ^ 1 order 
of magnitude more energy should be available in hadrons 
than electrons to produce these fluence ratios, when ac- 
count is taken of losses to neutrinos and inefficient ex- 
traction of proton energy. The neutral beam model inter- 
pretation therefore requires that the accelerated particles 
in, at least, GRB 941017 are hadronically dominated by 
^ 1 order of magnitude, consistent with predictions link- 
ing cosmic rays to GRB sources I Wick et al. (2003)] ). 

The neutron beam can take up to « 30 - 50% of the 
accelerated proton energy, with an equal amount of power 
going into neutrinos, gamma rays and electrons. At later 
stages when the internal photon field disappears, the cas- 
cade induced by electrons in the blast wave is suppressed 
as a consequence of the strong Klein-Nishina effects on 77 
photo-absorption for the external radiation, which is con- 
tributed mostly by backscattered prompt X-rays and MeV 
7 rays. But the backscattered radiation field still remains 
effective for photohadronic processes with UHE protons 
that continue to be accelerated in the GRB blast wave af- 
ter the decline of the prompt emission. Note that such 
an environment satisfies requirements for acceleration 



by external shocks ( Vietri, De Marco, fc Guetta (2003)] ), 
and for acceleration through the converter mechanism 



(Dcrishev, et al. (2003) I. Because the opacity to 77 ab- 
sorption is small for UHE 7 rays from ir° decay, these 7 
rays will escape from the blast wave and contribute to the 
neutral beam power at later stages directly. 

In conjunction with the neutral beam power is the 
prediction of a very significant fluence of UHE neu- 
trinos. From GRBs with anomalous 7-ray components 
and large fluences (the fluence of GRB 941017, includ- 
ing the anomalous component, is ^ 6.5 x 10~ 4 ergs 
cm -2 ), we predict several events detectable with 



In summary, we have proposed a hadronic model 
for the anomalous radiation component in GRB 941017. 
During the prompt phase, this emission originates from 
cascade radiation within the GRB blast wave. The anoma- 
lous component in the extended phase is synchrotron ra- 
diation from hyper-relativistic electrons formed by the as- 
sociated neutron beam escaping from the blast wave. The 
radiation components are well-suited for observations with 
the GBM and LAT on GLAST. GRBs with anomalous 7- 
ray emission components should also be bright neutrino 
sources detectable with IceCube. 
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